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A review of the Contractor Renormalization (CORE) metliod, as a systematic derivation of tlie 
ly-^ , low energy effective hamiltonian, is given, with emphasis on its differences and advantages over 

■ traditional perturbative (weak/strong links) real space RG. For the low energy physics of the 2D 

' Hubbard model, we derive the plaquette bosons (projected SO (5)) model which connects the mi- 

, croscopic model to phases and phenomenology of high-Tc cuprates. For the S = 1/2 Pyrochlore 

and Kagome antiferromagnets, the effective hamiltonians predict spin-disordered, lattice symmetry 
' breaking, ground states with a large density of low energy singlets as found by exact diagonalization 

\ of small clusters. 

l> 

CN ; I. RELIEF FROM STRONG FRUSTRATION 

' \ Frequently, interesting models of condensed matter systems involve strong local frustration. For example: the 
I ' Heisenberg antiferromagnet given by 

C/3 



h = jY,s,-s„ (1) 



I where {ij) are nearest neighbor bonds on lattices depicted in FigQ] The classical (infinite spin size) groundstates of 
the Pyrochlore and Kagome lattices, are known to exhibit macroscopic (exponential in system size) degeneracy, which 
^ can be lifted by quantum fluctuations. 

O \ At low enough temperatures, one expects the third law of thermodynamics to 'kick in' and that quantum fluctu- 
ations will choose a particular ground state. That ground state may, or may not, break spin rotational symmetry. 
However, sorting it out by semiclassical expansions such as spin wave theory, is a poorly controlled endeavor. In 
addition, numerical methods generally suffer from finite size limitations, and/or minus signs in quantum Monte Carlo 
^ ■ simulations. 

, Frustration causes fierce competition between nearly degenerate variational states, and equally plausible mean field 
' theories. Hence the phase diagram of such models are often a source of intense controversies. We advocate that 
\ such problems are best attacked by the 'divide and conquer' approach within a systematic real-space renormalization 
. scheme. The physical analogy is the use of nucleons, and then atoms, to treat the low energy correlations of the 
' standard model. It is obvious, that questions in chemistry, such as the relative stability of molecules, are better 
, resolved using effective interactions between atoms rather than by variational approximations on the high energy 
■ ('microscopic') interactions. 

, Returning to condensed matter physics, we have adopted the approach invented by Morningstar and Weinstein 
^ ' called Contractor Renormalization (CORE)P| to treat the square lattice Hubbard model 2||, and several problems of 
■i \ frustrated quantum antiferromamets|3, 0|- Other groups have also applied CORE to spin ladders [^, t-J ladders^, 
^ . and frustrated antiferromagnets 

Q ' The essence of CORE is that the microscopic lattice Hamiltonian is mapped onto an effective Hamiltonian which 
^ ,. acts on sites of a superlattice, within a lower energy Hilbert space, as we shall review below. After an effective 
^ hamiltonian is found numerically, and represented in terms of familiar second quantized operators (bosons, fermions, 
pseudospins etc.) the remaining task is to determine its ground state and excitations. This could be carried out in 
different ways. If the effective model is still highly frustrated (as we shall find for the pyrochlore case), the CORE 
H \ method could be reiterated. If the effective model turns out to be 'simple', that is to say: apparently unfrustrated, 
. - - . it naturally lends itself to variational approaches, and quantum Monte-Carlo methods (as was done for the Projected 
SO (5) theory of the square lattice Hubbard model|^, 0, and the Quantum Clock model of the Kagome 4J.) 

We shall see that if the effective interactions produced by CORE were calculated to all ranges, upto the size of 
the full lattice, the resulting effective hamiltonian would reproduce the exact spectrum of the original Hamiltonian. 
However, this in itself does not yield saving of numerical effort. The success of CORE relies on a rapid decay of 
the effective interactions with range. This range, which is derived from the numerical convergence tests, describes 
the physical coherence length of the effective degrees of freedom used to describe thee effective hamiltonian: e.g. 
bosonic hole pairs, for the square lattice Hubbard model, or pseudospins for the local singlets of the pyrochlore model. 
Therefore, a proper choice of effective degrees of freedom is helpful for rapid convergence. 
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FIG. 1: Strongly Frustrated Quantum Antiferromagnets treated by CORE. Red arrows denote classical spin directions. These 
can collectively rotate freely in the classical ground state manifold, rendering a poorly controlled spin wave expansion. 



II. CORE 



CORE is a non-perturbative block-spin rcnormalization, which uses exact diagonalizations to extract the effective 
interactions. 

1. Defining the reduced Hilbert Space. We first choose the elementary blocks which cover the lattice. (See Fig^for 
illustration for a square lattice). In order to preserve as much as possible the lattice symmetries of the original 
model (a choice of covering always breaks some translational symmetry), an optimal choice would be blocks 
which have the original rotational symmetries: such as plaquettes in a square lattice, triangles in the Kagome 
and triangular lattices, and teterahedra in the pyrochlore lattice. 

We diagonalize 7i on a single block and truncate all states above a chosen cutoff energy. This leaves us with the 
lowest AI states { | a) jf^. The reduced lattice Hilbert space is spanned by tensor products of retained block 
states I ai, . . . , ttAT ) . A case in point is the Hubbard model spectrum on a plaquette, which for the half filled 
case has 70 states (see Fig lHIII . We truncate 66 states and keep the ground state and lowest triplet, i.e. M = 4. 
Thus, the Hilbert space is considerably reduced at the first step. The retained states are in essence the 'atoms ' 
of the new effective Hamiltonian. The next task is to find their effective mass and interactions by calculating 
the intersite interactions. 

2. The Renormalized Hamiltonian of any cluster. The reduced Hilbert space on a given connected cluster of N 
blocks is of dimension M = . (See Fig. [21 for an illustration for N=3). We diagonalize Ti, on the cluster and 
obtain the lowest M eigenstates and energies: ( | n ) , e„), n = 1, . . . , A^. The wavefunctions | n ) are projected 
on the reduced Hilbert space and their components in the block basis | ai, . . . , ajq ) are obtained. The projected 
states ■i/'n are then Gramm-Schmidt orthonormalized, starting from the ground state upward. 

" \ m<n / 

where Z„ is the normalization. The renormalized Hamiltonian is defined as 

M 

^'■^"^^eJVin) (^J, (3) 

n 

which ensures that it reproduces the lowest AA eigenenergies exactly. 

Representing TY'"'^" in the real space block basis | ai, . . . , un ) , defines the (reducible) inter-block couplings and 
interactions. 
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FIG. 2: Covering the square lattice with plaquettes as elementary blocks. The reduced Hilbert space is defined as the tensor 
products of the lowest M states in each block. 



„ren 

(1,2,3) 



M 



states — ""^^ — ^"^^X^^ 




FIG. 3: A Cluster of 3 blocks, defining the effective Hamiltonian Ti.^"" of range 3. 



3. Cluster expansion. We define connected N point interactions as: 

(n, ■■■,«'«) 

where the sum is over connected subclusters of (ii, . . . , iat). The full lattice efi'ective Hamiltonian can be 
expanded as the sum 

i {ij) {ijk) 

hi is simply a reduced single block hamiltonian. hij contains nearest neighbor couplings and corrections to 
the on-site terms hi. hijk contains three site couplings and so on. hi-^^,,,^ij^ will henceforth be called range-N 
interaction. We expect on physical grounds that for a proper choice of a truncated basis, range-N interactions 
will decay rapidly with N. This expectation needs to be verified on a case by case basis. 
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In general, there is no a priori quantitative estimation of the truncation error. Nevertheless, if it decays rapidly 
with interaction range, we deduce that there is a short coherence length related to our local degrees of freedom, e.g. 
in our case the hole pair bosons and the triplets (bound states of two spinous). 

A. Comparison to perturbative real-space RG 

Perturbative real-space Renormalization of Quantum Many-Body systems is carried out in either the Lagrangian or 
Hamiltonian formulation. The Lagrangian renormlization involves integrating out of the path integral high wavevector 
and frequency modes V-'ingh- For example, in a model with point interactions of strength g, the renormalization is 
carried out by expanding the exponential in powers of g: models 

L^'^'li^io^] = i'V'highexp(-V*i(2)v>-5|Vp) 

« V*(a;,T)(L(2)+s(r-r'))v(a;',T')+/''"|Vl' + 0(/) (6) 

This formulation always truncates higher order terms in g (loops). It also necessarily introduces time-retarded 
interactions. This procedure usually results in a Lagrangian which similar to the microscopic one but with renormalized 
coupling constants. This allows an iterative renormalization (the renormalization group). However time retardation, 
if not neglected, divorces the Lagrangian formulation from the operator Hamiltonian formulation. 

The alternative is a perturbative Hamiltonian renormalization scheme. However, the traditional (non CORE) 
approach does not avoid the limitations of perturbation theory and time retardation. First one the Hamiltonian is 
separated into block terms (Hq) and inter-block interactions {H'). 

H = Ho + H' (7) 

The second step is to write the effective two-site hamiltonian in terms of a Brillouin-Wigner (BW) perturbation theory 
in the same reduced Hilbert space as CORE given by by the tensor products of the truncated block states \ai). 

where 

M 
i ai = l 

The expansion for Ti^^ contains intercluster interactions of all ranges, and the sizes of the terms is controlled by 
H' /A where A is a typical gap energy in the spectrum of Hq. The appearance of E inside Ti^^ (E) is a signature of 
time retarded interactions. It means that the spectrum is not given by the eigenvalues of H^^(£'o) for any choice of 

In summary, CORE has two major advantages over traditional perturbative real-space renormalization schemes: 

1. CORE is not an expansion in weak/strong bonds between block-spins. Its convergence does not necessarily 
depend on existence of a large gap to the discarded states of the Hilbert space. 

2. CORE is based on an exact mapping from the original Hamiltonian to an effective Hamiltonian, whose truncation 
error can be estimated numerically. Non-Hamiltonian retardation effects are avoided. 

III. SQUARE LATTICE HUBBARD MODEL 

An important interacting many-body model, especially in the context of high temperature superconductors, is the 
square lattice Hubbard model given by 

si 

^ = E ('^l^S-- + H-c) +uY^ ^n^ih (10) 



(8) 
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FIG. 4: The low spectrum of the four site Hubbard modeL Eigenstates are labeled by total spin S and plaquette momentum 
Qx, Qy = 0, TT. The undoped ground state is the vacuum, and the excitations are labeled by boson and fermion operators acting 
on that vacuum. 

where c|^, riis are electron creation and number operators at site i on the square lattice. Following the CORE procedure 
we choose to cover the square lattice by plaquettes (as in FigQ. The low spectrum of 2, 3, and 4 electrons (2, 1 and 
no holes respectively) is depicted in Fig lllll 

The ground state of the 4-site Hubbard model at half filling (rie — 4) is called | ) . which corresponds at large U/t 
to the resonating valence bonds (RVB) ground state of the Heisenberg model plus small contributions from doubly 
occupied sites. The product state 1 57) = Hf'"'' I ^) ji our vacuum state for the full lattice, upon which Fock states 
can be constructed using second quantized boson and fermion creation operators. 

The magnons are defined by the undoped plaquettes which are in the lowest triplet of 5* = 1 states. 

The hole pair state at {n,, — 2) is described by 

= ;^^i,o)Ti,o)ilO) 

= -ym E'^*j^»TCji + --- \^)^ (11) 



where dij is +1 (-1) on vertical (horizontal) bonds, and ... are higher order /7/t-dependent operators. Thus, 
creates a pair with internal d-wave symmetry with respect to the vacuum. For the relevant range of U/t, the state 
normalization is 1/3 < < 2/3. The important energy to note is the pair binding energy defined as 

Ab = E{0)+E{2)-2E{l) (12) 

where E{Nh) is the ground state of Nh holes. In the range U/t (0,5), it is bounded by — 0.04t < A;, < 0. It has 
been well appreciated that the Hubbard, t-J and even Cu02 models have pair binding in finite clusters starting with 
one plaquette. 
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A d-wave superconducting state can be written as the coherent state 

'plaq 

^d-scF ^ -Q ^^Qg Q ^ ee''fb^^) I n ) , (13) 

i 

with the superconductor order parameter 

mdijd^CjlY^) = ^Jz'^e'f sin 61 cos 6*. (14) 

Both the triplets and the hole pairs are 'bosonic states', which can be represented by boson creation operators acting 
on the RVB vacuum. They do not carry a negative sign under exchange. 

The single hole (3 electrons) ground states are fermions. Since they are slightly higher in energy than the hole 
pair states, we truncate the spectrum below them (at our peril, of course!), for the sake of deriving a purely bosonic 
effective hamiltonian, with hopefully rapidly decreasing interactions at long range. 



A. The Plaquette-Boson, (Projected SO (5)) Model 

We present the CORE calculations to range-2 boson interactions, while projecting out the fermion states. This 
required a modest numerical diagonalization effort of the Hubbard model on up to 8 site clusters. The resulting 
range-2 Plaquette Boson (PB) model can be separated into bilinear and quartic (interaction) terms: 

^rC[h] + n'[t] + w'[b,t] (15) 

where the bosons obey local hard core constraints 

h\h + Y,ti^U^<l (16) 

The bilinear energy terms are 

W*^ = (eb-2M)^6,K- Jfc5](6l6^.+H.c.) 

ia Oi{ij) 

-^E(^-4+H-c.)- (17) 

In Fig. IIII Al we compare the magnitudes of the magnon hoppings Jt , Ju and the hole pair hopping Jf, for a range 
of U/t. The region of intersection near U/t — 8, is close to the projected SO (5) symmetry point. We emphasize 
that although there is no quantum SO (5) symmetry in iJP*^, there is an approximate equality of the bosons hopping 
energy scales. This equality which was assumed in the pS0(5) theory |8||, previously appealed to phenomenological 
considerations. Here, the equality emerges in a physically interesting regime of the Hubbard model and has important 
consequences on the phase diagram as was shown in Ref . . 



B. CORE convergence and Coherence Length 

By diagonalizing the 12 site clusters, we have found that range 3 interactions are indeed between 1-10% of the 
range 2 interactions. Computing range 3 interactions /ii23 and finding out whether they are significantly smaller 
than range 2 terms is important for two main reasons: (i) This is the only way one could validate a truncation of 
the cluster expansion to range 2 for further investigations of the low energy properties of the model, and (ii) a rapid 
decrease in effective interactions signals a short 'coherence lengthscale' which describes the size of the effective degrees 
of freedom. For cuprates, the effective size of the hole pair is of experimental importance, since it is bounded by the 
superconducting coherence length as observed by the vortex core size, and the short superconducting healing length 
near grain boundaries and defects. Both have been observed to be not much larger than a few lattice constants. 
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FIG. 5: CORE results for the effective boson hopping rates corresponding to Hubbard interaction U (in units of t). Jb is the 
d-wave hole pair hopping rate, and Jt and Jtt are magnetic energy scales defined in 11711 . Note near equality of hopping energies 
around U/t = 8, signaling the projected SO(5) point. 



IV. PYROCHLORES 

The pyrochlore lattice is depicted in Fig^ Depicted in Fig. llVl is the spectrum of the Quantum Heisenberg model 
on a 4-site tetrahedron and a 16 site supertetrahedron. Both blocks can be used to cover the Pyrochlore lattice. 
Their block ground states are doubly degenerate singlets, and thus the effective hamiltonians are readily described by 
pseudospin-half operators. 

Using CORE for the tetrahedra covering upto range 3, we have found 

i?""" =E f(^(S--Sl)(S.-S)- (18) 

(Ijk) \ 




The coupling parameters (in units of J) are: J2 — 0.1049, J3 — 0.4215, and 0123, i = 1,2,3 are three unit vectors 
in the x-y plane whose angles 0^123 depend on the particular plane defined by the triangle of tetrahedral units 123 as 
given in table I of The effective hamiltonian 1118^ resembles the terms obtained by Tsunetsugu by second order 
perturbation theory (in inter-tetrahedra couplings) [13 ■ The classical mean field ground state is three of the four 
FCC sublattices are ordered in the directions e(0), e(27r/3), e(— 27r/3), while the direction of the fourth is completely 
degenerate. Therefore, classical mean field approximation for (|18|) is insufficient to remove the ground state degeneracy. 
Tsunetsugu^3 ^-ble to lift the degeneracy by including spinwave fluctuations effects which produce ordering at 
a new low energy scale. 

Here we avoid the a-priori symmetry breaking needed for semiclassical spinwave theory, by treating H18() fully 
quantum mechanically. This entails a second CORE transformation which involves choosing the "supertetrahedron" , 
as a basic cluster of four tetrahedra. 

Our new pseudospins are defined by the two degenerate singlet ground states of the supertetrahedron. (This 
degeneracy is found for the Heisenberg model on the original lattice as well as for the effective model (|18f) '). These 
states transform as the E irreducible representation of the tetrahedron (T^) symmetry group, similarly to the singlet 
ground states of a single tetrahedron. 

The supertetrahedra form a cubic lattice. The effective hamiltonian 118|l and the lattice geometry imply that 
non-trivial effective interactions appear only at the range of three supertetrahedra and higher. Range three effective 
interactions include two and three pseudospin interactions, which are dominated by 
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FIG. 6: The low spectrum of the Heisenberg model on a tetrahedron, and a supertetrahedron. In both clusters the ground 
states are degenerate singlets which can be represented by spin half eigenstates. 



^Cubic = J, ^ (r, . f,^. ) (r, . fy ) + (19) 

{(y)> 

Here, ( ) and (( )) indicate summation over nearest- and next nearest-neighbors, respectively. The couphng constants 
are found to be relatively small: Ji — 0.048J, = — 0.006J and J^'^ = 0.018J . The vectors fy depend on the 
vector Vij connecting the two sites, and their values are presented in Ref.^]. 

We performed classical Monte Carlo simulations using the classical (large spin) approximation to (|19|l . The ground 
state was found to choose an antiferromagnetic axis, and to be ferromagnetic in the planes as depicted in Fig. 
IIVI It differs from the semiclassical ground stateilOj. The latter involves condensation of high energy states of the 
supertetrahedron in the thermodynamic ground state. Since on a supertetrahedra we find a much larger gap to these 
states than intcr-sitc coupling, we believe they cannot condense to yield the semiclassical ground state symmetry 
breaking. 

To estimate the truncation error we calculated the contribution of range-4 interactions in both stages of CORE 
leading to l(TH|) and l(T^ . Evidently, these terms are small (< 30%), and most importantly, including them does not 
alter the mean field solution. 
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FIG. 7: Two CORE steps to relieve frustration in the Pyrochlore model. The first step results in an FCC pseudospin model 
whose mean field solution has one sublattice of completely free spins. The second CORE step results in a simple cubic model, 
with renormalized coupling constant of ~ O.OIJ, which has no residual macroscopic degeneracy. Its pseudospins, (which are 
represented by singlets on 16-site super-tetrahedra) , are antiferromagnetically correlated between neighboring planes. 




FIG. 8: CORE on the Kagome. Up-triangles provide four degenerate groundstates, which can represent a spin and a pseudospin 
of sizes half on each block. 



V. THE KAGOME 



The Heisenberg model (QJ on the Kagome lattice (depicted in FigP) has macroscopic degeneracy in its classical 
ground state. For the initial stage of CORE, we choose the upward triangles covering, and a truncated basis of the 
four degenerate spin half ground states, discarding the higher 5 = 3/2 states, see FiglVl. 

The S-L representation of the four ground states are labeled by |s, ^), where s =t, i is the magnetization and I =f|-, ^ 
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is the pseudospin in the z direction. Exphcitly, in the Ising basis IS1S2S3), 

{\^U)-\si\)) 



\sA) = 



V2 



I ||. 1-'^ Ti) + 1-^ it)) /2 

VA) = V 3'^^*'^*^' ^ ' 

The pseudospin direction in the xz plane correlates with the direction of the singlet bond, e.g. ff describes a singlet 
dimer on the bottom (— z) edge. Thus, the V eigenstates have definite chiralities. 



A. The SL Hamiltonian 



. The effective interactions between triangles is calculated by CORE. We note that this approach is feasible when 
two conditions are met: (i) Interaction matrix elements fall off rapidly with range such that the truncation error at 
finite ranges is small, and (ii) the norms of the projected eigenstates are sufficiently large for numerical accuracy. We 
have computed all range 2 and range 3 interactions, and neglected range 4 corrections, whose expectation values were 
found to be an order of magnitude smaller. At range 3, norms of projected eigenstates were greater than 0.75, with 
most states above 0.9. 

The effective Hamiltonian is a Spin-Pseudospin (SL) Model on the triangular lattice: 

'HsL = 'Hss + 'Hu 

(ij) 

+ Jssh(L,/" • Lj") + Jssiel(Li ■ Gij) 
+ Jssle2{^'j ■ Gji) + Jsslylylifiyj] 

+ JlylyVm (21) 

Here L-'^ = (L^,L^), and e^, e'^ are unit vectors in the xz plane. Tiss describes interactions of the Kugel-Khomskii 
type, where the pseudospin exchange anisotropy depends on the sites and bond directions. For any other other bond 
{if), Giji is simply found by rotating e^j by 0, ±120° according to the 0(2) rotation of {ij) {if)- 



B. Ground state and Low Excitations 



The best variational candidate for the SL model (|21|l are the dimer coverings of two triangle singlets, whose 
correlations are defined by Fig. IV Bl The dimer singlet states have been shown by Mila and MambrinijlJ to span 
much of the low singlet spectrum in finite cluster calculations. The variational analysis highlights the special role of the 
"Dimerization Fields", Jssiei, Jssiei in H21|) . for the formation of local singlets. These terms cancel under summation 
in all uniform states defined by (8^8^) = const. Their significant magnitude helps to lower the energy considerably by 
aligning with the anisotropy vectors eij to form singlets on certain bonds and not on others (Si • Sj) = — j^i^ij)^^^- 
This is a strong argument in favor of a paramagnetic ground state. Consequently, Tin is crucial in selecting the true 
ground state among the multitude of dimer singlet coverings. We have found that the perfectly ordered columnar 
dimer (CD) state minimizes Hu. A local "defect" of a rotated dimer in the CD background costs a "twist" energy of 
+0.01 per site. In (jj]) a theory of long wavelength fluctuations about the CD state has been investigated. The theory 
has a 6-fold 'clock mass term' ueCOs{64>), which yields a finite gap. However, quantum fluctuations renormalize down 
the magnitude of the clock mass to exponentially low values due to the 6 fold symmetry. This may explain the large 
density of low energy singlet excitations observed numerically and predicts a temperature dependence of the 
specific heat due to singlet excitations. 



VI. SUMMARY 



In summary I have reviewed our group's recent applications of CORE to currently interesting problems of strong 
quantum frustration, e.g.: the 2D Hubbard model for cuprates, and the geometrically frustrated Heisenberg model. 
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FIG. 9: Variational minimum of the SL Hamiltonian. Arrows denote pseudospin directions, and solid lines denote singlet spin 
correlations between neighboring blocks. The ground state is a columnar dimer state of neighboring up-triangle singlets. 

We find that the effective Hamiltonians are not necessarily simpler in form, but in a many cases they reveal the low 
energy degrees of freedom as the eigenstates of small clusters. The effective interactions, if they decay rapidly in 
space, may yield less competition and frustration than in the microscopic Hamiltonian, and thus be better amenable 
to variational solutions. 

In the computational sense, one could view CORE as an efficient algorithm to obtain the low energy physics of 
a large many body system, which maximally extracts its information from exact diagonalizations of small clusters. 
In combination with other approaches, this provides a promising direction to disentangle the low energy physics of 
strongly correlated many body problems. The method allows self-estimation of convergence, by sampling interactions 
of higher ranges than are retained. 
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